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Recently obtained field data and theoretical studies show that, for a fixed 
subscriber calling rate, dial-tone delay in No. 5 crossbar can be appre- 
ciably increased by high average line-link frame loads and also by poor 
load balance. The increased delay is caused by waste dial-tone-marker 
usage generated by a small number of calls that encounter horizontal 
group blocking in obtaining a dialing connection. This paper discusses 
an analytical model to predict the time from receiver off-hook to receipt of 
dial tone under various service conditions in No. 5 crossbar. 

I. INTRODUCTION 

1.1 Description of the dial-tone connection process 

In a No. 5 crossbar switching machine, dial tone is provided to a 
subscriber line, terminating on a line-link frame (llf), by an originat- 
ing register (or), terminating on a trunk-link frame (tlf), via a series 
of three network links : line links, junctors, and trunk links. The dialing 
connections are set up by dial-tone markers (dtms), which are common 
control devices. Each line-link frame contains a number of crossbar 
switches that are used to establish connections between subscriber 
lines and trunks, or between subscriber lines and service circuits, such 
as originating registers. The crossbar switches that form line concen- 
trators on which groups of subscriber lines terminate are called hori- 
zontal groups. Maximum size offices typically contain from 4 to 6 
dtms, 40 to 60 llfs, 20 to 30 tlfs, and up to 140 ors. Each llf contains 
10 horizontal groups and each horizontal group is a concentrator con- 
taining between 29 and 59 subscriber lines on the input side of the 
switch and 10 line links on the output side of the switch. 

To provide dial tone to a subscriber, the off-hook signal from the 
subscriber line initiates a bid for a dtm through a connector circuit 
unique to each llf. As soon as a dtm becomes available, it locates an 
unoccupied or and then attempts to find a dialing path (consisting of 
a line link, a junctor, and a trunk link) connecting the or with the 
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subscriber line. As soon as the connection is established, the dtm 
releases and proceeds to serve other calls waiting for dial tone. The 
or provides dial tone, receives the dialed digits, obtains a completing 
marker, transmits information to it, and releases. The completing 
marker then establishes the connection between the calling subscriber 
and an outgoing or intraomce trunk. 

If all ors are busy, the dtm releases and the call rejoins what is 
effectively a queue of calls waiting for a dtm. If an or is available, but 
no (unoccupied) dialing path connecting the or with the subscriber 
line can be found, the dtm will release the or, obtain a second (usually 
different) or, and try to find a dialing path between that or and the 
subscriber line. When a dtm is unable to find a dialing path between 
a given or and a given subscriber line, a matching failure is said to 
occur. If, on the second try, the dtm cannot find a dialing path, the 
dtm releases and the call rejoins the queue of calls waiting for a dtm. 
In such a case, a dtm second-failure-to-match (dtm2ftm) is said to 
occur. 

The method of assigning dtms to waiting calls is controlled by a type 
of "gating" circuitry designed to equalize service and to reduce the 
incidence of long delays in obtaining dial tone. When one or more 
dtms are free (light traffic operation), the llfs look for dtms according 
to a fixed preference order. When all dtms become busy and a request 
for a dtm occurs, a gate is closed and only the llfs containing calls 
waiting for dial tone at that moment are put inside the gate. The 
dtms then proceed to serve the llfs inside the gate. If more than one 
call is waiting on an llf, only one call will be served during the gating 
period. Once a call on an llf is served during a given gating period, 
that llf is put out of the gate, whether or not the dtm is successful 
in establishing a dialing connection for the call. When all llfs with 
requests at the start of the gating cycle have been put out of the gate, 
the gate opens ; if there are sufficient waiting calls to cause all of the 
dtms to become busy again, a new gating cycle will be started ; other- 
wise, light traffic operation will resume. During a gating period, a dtm 
that becomes idle scans the llfs in cyclic order. Each dtm has a 
different starting llf for the scan so as to equalize service. (The de- 
scription of the gating procedure is taken from Refs. 1 and 2.) 

1.2 Effects of matching failures on dial-tone delay 

Repetitive matching failures can occur in establishing the dialing 
connection essentially because the first-stage crossbar switch on which 
the subscriber line terminates (the horizontal group) is a concentrator 
whose output links (line links) have holding times that are much 
longer than the holding times of the common control devices that set 

1756 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1975 



up the dialing connections (the dtms). The average line-link holding 
times, being largely determined by conversation holding times (with 
allowance for ineffective attempts), can be on the order of 150 seconds 
or more, while the dtm holding times are typically on the order of 
0.25-0.40 second. 

Although matching failures (blocking) can occur when one or more 
of the 10 line links on a horizontal group are unoccupied, repetitive 
matching failures under such conditions are rare because a dtm is 
quite likely to be successful in establishing a dialing connection after 
a few attempts. The expected holding time in a blocked condition of 
a call that finds all 10 line links busy is the lesser of the length of time 
for one of the 10 line links to become free (with a short added time 
for dtm uses and matching failures in setting up the dialing path) and 
the length of time that a subscriber is willing to wait. 

Since, at average busy-hour load levels, less than about one percent 
of all originating calls are predicted to encounter an all-10-line-links- 
busy condition, previous dial-tone-delay studies have assumed the 
effect of matching failures on dial-tone delay to be small. However, 
the fact that line-link holding times are much longer than dtm holding 
times means that a call which finds all 10 line links busy can remain 
blocked for long enough to consume a large number of dtm uses (except 
at calling rates sufficiently high that very little of the offered blocked- 
call load is carried by the dtms). Thus, it is possible for a small number 
of calls experiencing blocking to generate a disproportionate number of 
waste dtm uses, increase dtm occupancy, and thereby increase dial- 
tone delay for all other calls in the office. 

II. ANALYTICAL MODEL (LIMITING FORM) 
2.1 Assumptions 

Since the main effect of matching failures on dial-tone delay is 
caused by the resulting waste dtm use, and since the dial-tone delay 
distribution of the small proportion of calls experiencing repetitive 
matching failures can be calculated approximately (see Section 5.2), 
the effect of matching failures on the dial-tone delay distribution of 
calls that do not experience repetitive matching failures has been repre- 
sented in terms of a queuing model with two classes of calls, good calls 
and bad calls, defined as follows : 

(i) A good call experiences no matching failures but is subject to 
delay caused by dtms and ors. 

(it) A bad call experiences total network blocking (no dialing con- 
nection available) and defects from the system after an expo- 
nential waiting time. 
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We will first describe the mathematical structure of the analytical 
model used for dial-tone delay calculations. We will refer to this model 
as the limiting model. Next, we derive the dtm saturation load* and 
prove that it is not changed by the presence of bad calls. 

In Appendix B, we prove that the equilibrium queue length and wait- 
ing time distributions of the limiting model are the limits in distribu- 
tion of a sequence of equilibrium distributions arising from a model in 
which the expected bad-call arrival rate approaches zero and the 
expected bad-call waiting time (until defection) approaches infinity 
in such a way that their product, total erlangs of bad calls, is constant. 

In the limiting model, a good call arrives and finds a random 
(truncated) Poisson-distributed (but time-independent) number of 
bad calls permanently present in the system. 1 " The queue discipline 
is characterized by random order of service. The dtms cannot dis- 
tinguish between good and bad calls when choosing a call to be served. 
This corresponds to the fact that, in the No. 5 crossbar switching 
machine, a dtm cannot recognize that a 10-line-links-busy condition 
exists on a particular horizontal group and, hence, cannot avoid wasting 
time serving calls for which no dialing path exists. Equilibrium good- 
call queue length distributions are computed conditional upon the 
number of bad calls present in the system. Calculating the expectation 
of the conditional distribution over the distribution of bad calls gives 
the unconditional dial-tone-delay distribution for good calls. The con- 
ditional distributions depend on the total office calling rate and on the 
number of llfs, dtms, and ors but do not depend on the horizontal 
group load or the load variation. Hence, delay distributions for a range 
of frame load and balance effects can be calculated using the same set 
of conditional distributions, thereby greatly reducing the computer 
time needed for parametric studies. 

Since dtm holding times are approximately constant (for a given 
set of office parameters and traffic characteristics) and since these 
holding times are small with respect to the accuracy with which it is 
necessary to be able to predict delays, we treat the dial-tone delay 
process as a discrete time queue with a constant service time of T 
seconds. Good calls are assumed to arrive in batches according to a 
Poisson process at times kT, for k — 1, 2, • • •. Immediately upon 
arrival, the good calls join the queue of good and bad calls waiting 
for dial tone. Calls are chosen at random from the queue for service 



* The dtm saturation load is defined to be the good-call originating load below 
which a steady state good-call-queue-length distribution exists and above which 
such a steady state distribution does not exist. 

+ When the bad-call input is assumed to be peaked, the number of bad calls in 
the system has a (truncated) negative binomial distribution. This is discussed in 
Section IV. 
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by the dtms with each call — whether a good call or a bad call — having 
an equal probability of being chosen. A good call served on the dtm 
cycle beginning at time kT will either acquire an or, receive dial tone, 
and thereupon exit the system at time (k + 1)T, or else will fail to 
acquire an or and will return to the queue of calls waiting for dial 
tone at time (k -f l)T. A bad call served on the dtm cycle beginning 
at time kT will return to the queue of calls waiting for dial tone at 
time (k + 1)T. Note that since new arrivals occur only at the time 
points kT, a dtm that is idle on the cycle beginning at time kT will 
remain idle at least until the start of the cycle beginning at time 
(k + 1)T. The dial-tone delay for a call is the length of time from the 
moment when the call arrives to when it obtains dial tone and (simul- 
taneously) leaves the dial-tone queue. Thus, the minimum possible 
dial-tone delay in the model is T seconds. 

In the analytical model, two dtms serving N line-link frames are 
used to represent four dtms serving 2N line-link frames. A single dtm 
holding time equal to the office average dtm holding time is used for 
both good and bad calls. In actuality and in the simulation models, 
bad calls have somewhat longer dtm holding times than good calls 
and good calls that encounter a condition of all-ORs-busy have some- 
what shorter dtm holding times than good calls that do not encounter 
a condition of all-ORs-busy. Comparison of results from the analytical 
model with those of the gating simulation model indicates that these 
simplifications tend to offset each other in the range of interest. 

A further simplification in the analytical model concerns the manner 
in which availability of ors is treated. Since No. 5 crossbar offices 
typically contain over 100 ors and since the ors are in tandem with 
the dtms, it is presently not practical to keep track of the number of 
occupied ors directly in an analytical model. Availability of ors is 
treated by assuming that at each time point kT, all ors are busy with 
probability q and two or more ors are free with probability p = 1 — q. 
The calculated probability that exactly one or is free and that a 
dialing connection is available between this one or and the given sub- 
scriber line is sufficiently small, in the occupancy range of interest, 
so as not to warrant the additional complexity caused by introducing 
this effect. 

The idea of using a discrete time model and the method of treating 
or availability through use of a fixed probability of all ors busy are 
due to Halfin. 3 Following Halfin, the probability q is taken to be the 
erlang C probability of all ors busy. At or occupancies above about 
0.90 with frame loads low enough that few second failures to match 
occur, this method of treating or availability somewhat underpredicts 
the delay caused by an all-registers-busy condition (based on com- 
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parison with simulation results). The underprediction arises because, 
at higher or occupancies, once all ors become busy they tend to remain 
busy for a time period equal to several dtm cycles, as one would expect 
from the erlang C delay formula. 

Preliminary studies were made using an analytical single-server 
cyclic queuing model, developed by S. Halfin, 3 which represented the 
No. 5 crossbar gating process in considerable detail but did not take 
into account the effects of horizontal group blocking. These studies 
showed that the delays predicted by the cyclic queuing model do not 
differ appreciably from those of a discrete-time M/D/l queue with 
feedback and random order of service. The latter model requires 
about l/30th of the computer time required by the former. Both 
models overpredicted simulations. For these reasons, explicit repre- 
sentation of the gating process was not attempted and a discrete- 
time M/D/2 queue with feedback and random order of service was 
taken as the starting point for developing an analytical dial-tone-delay 
model to include effects of horizontal group blocking. 
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Fig. 1 — Comparison of analytical and simulation models. Dial-tone delays are at 
1400 ccs/llf and 0.70 dtm occ (excluding dtm2ftms). 
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Comparison of simulation results with predictions using the limiting 
analytical model (based on two markers and random order of service) 
shows that the predicted delays typically fall somewhat above or very 
close to simulated delays based on four markers with gating order of 
service and somewhat below simulated delays based on four markers 
with random order of service. In the light of these results and because 
of the large scatter in actual measured No. 5 crossbar dial-tone delays, 
it did not seem worth the considerable added complexity to include 
explicit representation of the gating process in the analytical model. 
Typical results are shown in Figs. 1 and 2. 

2.2 Queue length equations for the limiting model 

The queue-length process for good calls in the limiting model is a 
discrete-time Markov chain with finitely many irreducible classes 
{Cfc} that are noncommunicating in the sense that no transition be- 
tween different classes is possible. Hence, each class is itself an irre- 
ducible Markov chain. The fcth class is the queue-length Markov chain 
for a discrete-time M/D/2 queue with random order of service, 
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Fig. 2 — Dial-tone delays are at 1600 ccs/llf and 0.80 dtm occ (excluding dtm2ftms). 
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"feedback" (occurring when an all-ORs-busy condition is encountered), 
and k blocked calls residing permanently in the system. The dtms 
cannot distinguish these k permanently blocked bad calls from the 
good calls in the system, so each of the k bad calls and each of the good 
calls in the system compete on an equal basis for dtms. Upon comple- 
tion of service by a dtm, a bad call rejoins the queue. With probability 
q, the good calls served at any given time period rejoin the queue and 
with probability p they leave it. We will show that each of these 
queuing systems has the same saturation load and we will describe 
how the equilibrium-queue-length distribution and the equilibrium- 
waiting-time distribution for each system is calculated. 

Viewing the whole process again from the standpoint of a Markov 
chain with finitely many noncommunicating classes, the queue-length 
and waiting-time distributions referred to above may be regarded as 
being conditional on k, the number of blocked calls (permanently) 
present in the system. Now let k be a random variable with a truncated 
Poisson distribution of mean x, or, equivalently, regard the Markov 
chain as having any initial distribution ir(k, i), where i denotes the 
number of good calls in the system and where the marginal distribu- 
tion ir(k, •) of the number of bad calls in the system is truncated 
Poisson with mean x. Then the (unconditional) equilibrium-queue- 
length and equilibrium-waiting-time distributions for this system may 
be computed by taking the expectation (with respect to the truncated 
Poisson distribution of k) of the conditional queue-length and waiting- 
time distributions for each of the individual systems represented by the 
classes (7*. 

Let X„ denote the number of good calls in the queue at time nT and 
let Y denote the number of bad calls (permanently) present in the 
queuing system. Let 

P k (i, j) = Pr [X n+1 = j\X n = i and Y = A-] 

and let Pk(i) be the equilibrium-queue-length distribution for the case 
of k bad calls. Let A (n) be the probability that n good calls arrive in 
one service time interval (of length T) . Then, by assumption, 

A(n) = e ^f n- 0,1,2,..- 

= otherwise, 

where 

A = XiNT 

X x = llf originating call rate (calls per second on one llf), 
N = number of llfs served by the two dtms. 
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The transition functions P k ( ■ , • ) are given by 

Po(0,j) =A(j) 

Po(l, 3) = pA(j) + qA(j - 1) 

Po(i, j) = pA (j -i + 2)+qA (j - for i ^ 2, (1) 

and for k ^ 1, 
P*«W) =A(i) 

P*& i) 
= p[ft(*-iU(j-»)+2»fcA(i-t+i)+i(t-i)4(i-t+2)]/ 

[(* + *)(* + *- l)] + qA(j-i) for » £ 1. (2) 

Note that P k (i, j) = f or i > j + 2 and that P t (j + 2, j) > 0. 

If the equilibrium distributions P k (-) exist, then they satisfy the 
equation 

P*(j) = L Pk(i)P k (i, j). (3) 

Let f k be the generating function for the fclh queue length process ; 
then 

/*(*) = f Pk(i)z\ 

t— o 

The generating function / is easily seen to be given by 



f (a = p(* - DEC* + 1)^0(0) + ggoO)] 

• /0V ' AXO-.) - (p + ? 2 2) 

If P (-) exists, then lim z ti/ (2) = 1, and hence 



(4) 



Po(0) + ip (i) = i-A. (5) 

Let S k 2) (p) denote the saturation load of the fcth process with two 
dtms, as discussed above, and let Sjp(p) denote the saturation load 
of the analogous one-DTM system with k blocked calls. It can be shown 
that SIP(p) = p for all k = 0, 1, 2, • • •, i.e., a necessary and sufficient 
condition that a steady-state queue length distribution exists in the 
one-DTM system is that X < p.* Using this result, we will show that 
S k 2) (p) = 2p for all k = 0, 1, •••. 

Since the Markov chain for the kth process above is irreducible, 
P (0) > and P (l) > whenever the stationary distribution P (-) 
exists. Hence, it follows from eq. (5) that SJ 2) (p) ^ 2p. Since S k 2) (p) 



* This result can be proved using a theorem of Kushner (Ref. 4) and some properties 
of generating functions. 
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^ Sif(p) ^ S^(p) ^ 2p for k 2 ^ fcx ^ 0, it suffices to show that 
$ 2) (p) ^ 2pforallfc ^ 0. 

Consider a two-DTM system with 2k blocked calls, which differs from 
the 2fcth process defined above only in that (i) each dtm maintains a 
separate queue with half the originating traffic being assigned to one 
of the dtms and the other half being assigned to the other dtm and 
(u) each dtm serves k blocked calls. Let S 2 k(p) denote the saturation 
load of this modified system. Then clearly, S 2k (p) ^ 8jS(p) and, by 
symmetry, Snip) = 2SfP(p) = 2p, because the modified system 
simply consists of two identical one-marker systems working sep- 
arately, with each being assigned half the incoming traffic. Hence, 
8S(p) ^ 2p and so Si 2) (p) = 2p for all k ^ 0. 

Thus, for X/2p < 1, a stationary distribution exists for the fcth 
process for each k ^ 0. In the sequel, we will confine ourselves to the 
case where X/2p < 1. In this case, /o exists and is analytic for \z\ < 1 
and continuous for \z\ = 1, so the numerator of the expression on the 
right-hand side of eq. (4) must vanish for any \z\ ^ 1 for which the 
denominator vanishes. It is easy to see that the denominator has a 
single real root on the open interval ( — 1, 0) and that z = 1 is a root. 
By applying Rouche's theorem to the functions z 2 exp[\(l — z)~\ 
and p + qz 2 along the circle | z \ = 1 + e, one can show that for X/2p 
< 1 and for e > sufficiently small (depending on X and p), the ex- 
pression z 2 exp [X(l — z)~\— (p + qz 2 ) = has exactly two roots in 
the open disc \z\ < 1 + e. Thus, the root on ( — 1, 0) and the root 
at 1 are the only roots within the closed unit disc and 

(« + l)Po(0) + *Po(l) = 0, (6) 

where £ is the unique root of the expression 

2 2gX(l-z) _ ( p _|_ qz 2) = 0> -1< 2 <0. (7) 

Solving eqs. (5) and (6) for P (0) and P (l) yields the unique 
solution 

_ -€(2p-\) 
Po(0)_ (l-Qp 

Po(1) ~ (1 - Qp ' (8) 

where £ is denned by (7). The values P (j) for j ' £ 2 can now be com- 
puted recursively from (3) using (8) . 

When k > 0, computation of P*(-) by the method of generating 
functions involves numerical solution of some rather unwieldy 
differential equations. The following approach, which makes use 
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of properties of recurrent Markov chains, is easier to implement 
computationally. 

In Appendix A we show that, for X/2p < 1, P k (-) is given by 
c(o + r§), where c and r are constants and a, (3 are the two (particular) 
eigenvectors, defined by (9) below, of the transition matrix P k (i, j).* 
We further show that the condition, a, + rjSy ^ for each j, deter- 
mines r uniquely and provides an easy way to calculate r. Once r is 
obtained, the constant c is uniquely determined by the condition that 
Hf=aPk(j) = 1. Since the quantities ay and /3y can be computed re- 
cursively, the above observations lead to a simple algorithm for com- 
puting Pk(-). The results are given below, the proof is in Appendix A. 
Let 

OLQ = 1 

a, = 

n+l 

a» - E <*jPk(j, n) 
a " +2 = ~ Pt(n + 2,n) W = ° < 9a) 

/3„ = 
ft = 1 

ft. - E 1 ftP*0", n) 

^ = "- P^+2^)- - * °- ( 9b ) 

It is shown in Appendix A that there is a unique constant r satisfying 

m„ £ r ^ M n all n ^ 3 (10a) 

lim ra„ = r = lim M„, (10b) 

n-»co n-»oo 

where, for ?i ^ 3, the increasing sequence m n and the decreasing se- 
quence M n are defined by 



!„ = max - !*: a, < 0, ft >0 

f.--j{«[-g^>*A<o]- 



I .-• . ■ ■ - j ■ i 

(11) 



Since the quantities m n and il/„ can be computed recursively, eqs. 
(10a), (10b), and (11) yield a well-defined algorithm for computing r. 

Once r has been computed, P k (j) can be calculated from (3), setting 
P*(0) = c, P*(l) = re and determining c by the requirement that 

E?-o Pk(j) = l. 



* Each of c, r, a, and (J depends on t. For each integer j ^ 0, ay and /Sy, denote the 
jth components of the eigenvectors a and (3, respectively. 
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Since recursive calculation of Pk(j) by (3) involves successive sub- 
tractions, it is essential to perform all computations in double precision. 

In a No. 5 crossbar switching machine, the measured dtm occupancy 
during an hour is defined to be the total work time (in seconds) of all 
the dtms divided by the product of the number of dtms and the number 
of seconds in an hour. Thus, the measured occupancy is the average 
fraction of time that a given dtm is occupied. In the model the (steady- 
state) dtm occupancy is denned to be the limit, as n approaches in- 
finity, of the expected fraction of time that a given dtm is occupied 
during the n work cycles of length T beginning at times T, 2T, • • • , nT. 
This limit is equal to the probability that a given dtm is occupied on 
a work cycle beginning with the system in steady state. Hence, 

dtm occupancy = P (two or more calls are waiting for dial tone) 

+ \ P (exactly one call is waiting for dial tone) 
= 1 — P (no calls are waiting for dial tone) 

— \P (exactly one call is waiting for dial tone). 

As discussed earlier, the number of bad (blocked) calls waiting for 
dial tone at any given time has a truncated Poisson distribution in the 
model. Let 

x = total erlangs of blocked calls on two dtms 

K = maximum possible number of blocked calls that can be in the 
system (waiting for dial tone) at any time* 

XV /y. K 

c x k = L rv 

4 = K ! 

Then the probability that fc blocked calls are in the system is given 
by (x k /k \)cJk- Let p„ denote the (steady state) dtm occupancy for 
the case of no blocked calls in the system and let pb denote the (steady 
state) dtm occupancy including the effect of blocked calls. It follows 
from (5) that 

Pu = X/2p. (12a) 

Hpticp 

Pb = 1 - c7APo(0) - § [cJiPo(l) + xc£Pi(0)] 

= 1 - ci(l - Pu) - (x/2)c7&Pi(0). (12b) 

2.3 Delay equations for the limiting model 

In this section, we calculate the dial-tone delay probabilities for a 
call that arrives when the system (just prior to the arrival of the call) 



*The truncation parameter K determines the number of conditional delay dis- 
tributions to be calculated, so K should be taken to be no larger than needed to retain 
sufficient accuracy in the calculations. For values of x in the range of interest e~ x c x g 
differs from 1 by less than about 10 -3 for K = 5 ; thus, for computational purposes, 
K may be taken to be 5. 
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is in steady state. Recall that the queuing model is in discrete time 
with arrivals at times nT and dtms operating at times nT. Calls served 
on the dtm cycle beginning at nT either exit the system at time 
(n + l)T or else return to the queue at time (n + l)T. 

Thus, the steady-state probability that a call experiences a dial- 
tone delay of (m + 1)T seconds can be computed recursively in terms 
of the probability of a delay of mT seconds by straightforward condi- 
tional probability calculations involving matrix multiplications. 

Let 

Y = number of bad calls (permanently) present in the system 
X n = number of good calls in the queue at time nT 
X = number of good calls in the queue immediately after an arrival 
when, just prior to the arrival, the system was in steady state. 

Thus, X is the total queue length just after arrival of a call when the 
system is in steady state. Let 



Pk(j) = P 



Then 



X n = j 



Queue is in steady state at time (n — 1)T^\ 
with k bad calls permanently present in 
the system, and at least one good call 
arrives at time nT 



P[± = j\Y = k-] = P k (j). 



Let 



P k (i, j) = p (x n+1 = j 



X n = i, Y = k, and at least one good' 
call arrives at time nT 



and 



A(n) = 



A{n) 



1-A(0) 
= 



n ^ 1 

otherwise. 



Then A(n) is the conditional probability that n calls arrive at time 
nT given that at least one call arrives at time nT and Pk(i, j) is defined 
by eqs. (1) and (2) with A used in place of A. Also, 



y+2 



PUj) = E Pk(i)Pk(i, j). 



Let W n {i, k) be the conditional probability that a call arriving in 
steady state has a dial-tone delay of nT seconds, given that the queue 
length of good calls upon arrival of the call is i and that the number 
of bad calls permanently present in the system is k. Thus, letting D 
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denote the dial-tone delay, we have 

W n (i, k) = P[D = nT\X = i and Y = h~] i ^ 1. 
Then 

Wi(i, k) = min |^1, (fe + g J p » ^ 1, 

and PT„ + i(-,/e) can be computed recursively from W„(-,k) as de- 
scribed below. 

Let X n = i ^ 1 and consider any one given call out of the i calls 
present at time nT. Let 

Z B+ i = j and the given call does 
not leave the system on the dtm 
cycle beginning at time nT 



<P*tt 3) = P 
Then 



X„ = i and F = k 



W„ + i(i,k) = £ <P 4 & j)W n (j,k), 



and, for t «£ 1, (P*(*, i) is given by* 

<Po(», j) = [l - mln f ? ° ][M(j - t + 2) + flA(j - 0] 

+ [=*f^]«A<;-<> 



fl»i(», i) = 



t — 



i + 



■j ][?[(» " 2)A(i - i + 2) + 2A(j - i + 1)] 

+qA(j -**)] + [ i-|^ ] ff^W ~ 
and, f or k ^ 2, 

^ tt /> - [i - jht-J 






p £, >'+'<" ^ "-'°-' + ' ) 



4-gA(i - i) 



+ 



(iTl)^tf-«. 



The second term on the right-hand side of the above equation is 
the probability of the event that: (i) the given call is selected for 



* The derivation of <P k (i, j) for k ^ 2 is given below. The derivations of (P*(t, j) 
for k = and k = 1 are analogous and are omitted. 
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service by the dtms, (m) all ors are busy, and (Hi) X n+l = j. (Recall 
that, in the model, either all ors are busy, with probability g, or else 
two or more ors are available, with probability p = 1 — q). The 
first term is the probability that : (i) the given call is not selected for 
service and (ii) X n+ i = j. In the first term, the expression in large 
brackets is the conditional probability that X n +i = j given that the 
call is not selected for service. This conditional probability is itself 
composed of the probabilities of two mutually exclusive events. The 
second term within the large brackets is the (conditional) probability 
of the event that all ors are busy and X n+1 = j. The first term within 
the large brackets is the (conditional) probability that two or more 
ors are free and X n+1 = j. The sum from I = to I = 2 in this term 
pertains to the cases where 0, 1, or 2 good calls, respectively, are selected 
for service by the dtms. Since this sum is part of a probability condi- 
tional upon a given good call not having been selected by the dtms, 
the available population from which the dtms may select calls con- 
sists of k bad calls (ft j> 2) and i - 1 good calls (i ^ 1). Since the 
selection is without replacement, the selection probabilities have the 
hypergeometric form shown above. 
By the law of total probabilities, 

P£D m nT\ Y = ft] - f W n {i, k)PlX = i\Y = k], 
P[Z> > ?iT\ Y = kl = 1 - jr P[D = mT\ Y = k\ 



and 

prn -> nT\ = 

*tb ft! 



P[D > nT~\ = t C -4r P£D >nT\Y= ft} (13) 



where 



c xK = £ x k /k !. 

k=0 



III. CALCULATION OF OFFERED BLOCKED-CALL LOAD 

In the limiting analytical model, the number of blocked calls in the 
system has a time-independent truncated Poisson distribution with 
mean x. This section describes a method for computing x, the mean 
offered erlangs of bad calls, in terms of the distribution of carried load 
among horizontal groups in the office. Using the limiting analytical 
model in conjunction with these methods for calculating the offered 
blocked-call load, we can calculate the No. 5 crossbar dial-tone delay 
distribution in terms of the calling rate and the distribution of carried 
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load among horizontal groups in the office. Thus, we can predict the 
effect of frame load and balance on dial-tone delay in No. 5 crossbar 
offices. 

We first construct a model for an individual horizontal group and 
express the expected bad-call load from one horizontal group as a 
function of the carried horizontal group load. To get the total ex- 
pected bad-call offered load for the office, the expected contribution 
from an individual horizontal group is integrated (numerically) over 
the office horizontal group load distribution. The model, described 
below, for representing an individual horizontal group may be called 
a modified finite source Palm delay model. 

We assume that the input to a horizontal group is from a finite 
number of sources with equal calling rates and that call-holding times 
are exponential. The number of sources N is taken to be 35 although 
in actuality most horizontal groups have 49 or 59 subscriber fines. 
The reason for the use of the lower number of sources is that an earlier 
study by W. S. Hayward, Jr. 5 showed that blocking on concentrators 
with unequal line occupancies can be approximated by blocking calcu- 
lations based on equal calling rates and a lower number of sources. 
The use of 35 sources was suggested by J. G. Kappel. 8 

The calculations take into account the fact that an incoming call 
cannot occupy a subscriber line when all ten line links are occupied. 
Calls that find all line links busy will either defect or will eventually 
obtain a line link. While waiting for a line link to become available, 
a call is assumed to have an exponential waiting time until defection, 
with a mean of 30 s.* 

In the case of ideal load balance, each horizontal group in the office 
is assumed to have a true carried load of z erlangs. In the case of less 
than ideal load balance, the distribution of true carried load among the 
individual horizontal groups is assumed to be normal with mean z 
and coefficient of variation c g , where z is the office average carried 
horizontal group load and where c B is the group-to-group coefficient 
of load variation for the office. The term c„ may be inferred from office 
load balance data either by using analysis of variance or, more com- 
monly, by subtracting a standard value of the residual variance from 
the total measured variance of the office horizontal group load 
distribution. 



* This is the same value that was used in step-by-step dial-tone-delay calculations. 
(See Ref. 7.) These calculations are based on the Palm delay model 8 using an assumed 
mean call holding time of 150 s and an assumed j factor of 5. In the notation of Ref. 
9, this value of j corresponds to a mean-time-to-defection of 30 s. In Ref. 9, it is also 
stated that this value was found to be slightly conservative for most applications, 
based on review of panel office data reported in Ref. 9 and other (unpublished) 
step-by-step data. 
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To obtain the total offered blocked-call load for the office, we write 
the offered blocked-call load for an individual horizontal group x(z) 
as a function of the individual horizontal group carried load z and 
integrate the function over the office distribution of carried load. We 
now describe how x(z) is calculated. 

Let N denote the number of subscriber lines per horizontal group. 
As discussed above, N has been taken to be 35 in all computations. 
Let X n denote the combined originating and terminating rate when 
n subscriber lines are occupied and let n„ denote the subscriber line 
hang-up rate when n subscriber lines are occupied. Then 



X„ = X(iV - n) ^ n ^ 9 

X 
2 



= 7T (N - n) 10 ^ n ^ N - 1 



and 

Hn = n/H ^ n g 10 

= 10/H + (n - 10)/H b 11 S n £ N, 

where H denotes the office average line-link holding time and H b de- 
notes the reciprocal of the defection rate for a call that is waiting for 
a line link to become available. As discussed above, H b is taken to be 
30 s based on results in Ref. 9. 

The parameter X (combined originating and terminating rate per 
unoccupied subscriber line) is an unknown whose value will be ob- 
tained from the horizontal group-carried load z. The factor of f appear- 
ing in the definition of X„ for 10 ^ n ^ N - 1 reflects the fact that, 
when all 10 line links are busy, an incoming call cannot cause a sub- 
scriber line to be occupied. The definition of the hang-up rate /m for 
11 ^ n g> N reflects the assumption that, when all 10 line links are 
occupied, the holding times of the subscriber lines for which no line 
links are available should be shorter than full call holding times. 

For a horizontal group with carried load z, let 

7T„ = steady- state probability that n (out of A*) subscriber 
lines are occupied. 

Then, using standard methods for computing the steady-state dis- 
tribution of a birth-and-death process, 10 

iro = c and 7r„ = c ( ° " " -1 ) for n > 1. 

where the constant c is determined by the requirement that 

N 

X) T n = 1. 

n=0 
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Then the carried horizontal group load z is given by 

2 = £ n-TTn + 10 £ x B . (14) 

n=l n-10 

It is not hard to show that the carried horizontal group load z is a 
strictly increasing function of the subscriber line occupancy rate X and 
that z ^ min [10, \NH], This makes it easy to determine numerically 
the unique value of X corresponding to a given carried load z. 

Once X has been determined, the quantities X„, n„, and w n are used 
to compute the offered blocked-call load contributed by the hori- 
zontal group. We take x{z) to be the expected number of occupied 
subscriber lines for which no line links are available. This yields the 
value 

x{z) = £ (n - 10)x„. (15) 



IV. PEAKEDNESS OF THE BLOCKED-CALL STREAM 

In the model discussed in Section III, blocked calls arrive according 
to a Poisson process and defect after an exponential waiting time. 
Since blocked calls constitute an overflow stream and since it is well 
known 11 that overflow traffic usually has a peakedness* greater than 1 
and hence is not Poisson, some discussion of the peakedness of the 
blocked-call stream is in order. 

The blocked-call stream is the superposition of overflow traffic from 
all of the (typically 400 to 600) horizontal groups in an office. In the 
case of an office with ideal load balance (identical horizontal group 
loads), all horizontal groups would have equal expected contributions 
to the blocked-call stream and standard limit theorems would suggest 
that the blocked-call stream should be approximately Poisson. 

Comparisons of calculated and observed dial-tone delays discussed 
in Section V (covering measured dtm occupancies up to about 0.84) 
indicate that, when the blocked-call stream is assumed to be Poisson, 
the calculated delays generally fall in the midrange of applicable data. 
However, there is a large variability in observed dial-tone delays 
measured in the same office under nearly identical levels of dtm 
occupancy, second-failures-to-match, and percent all-ORs-busy. The 
presence of this variability may be regarded as evidence that, in some 
busy hours, the blocked-call stream may be peaked in nature. Peaked- 
ness of the blocked-call stream is capable of accounting for substanti- 



* The peakedness of a stream of calls is, by definition, the variance-to-mean ratio 
of the number of busy servers when the stream is offered to an infinite group of 
servers with independent identically distributed exponential holding times. 
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ally higher calculated dial-tone delay at a given level of dtm occupancy, 
second-failures-to-match, and percent all-ORs-busy than would be 
calculated under the assumption that the blocked-call stream is 
Poisson. 

One would expect the blocked-call stream to be peaked whenever 
most of the blocked calls are contributed by a small number of highly 
overloaded horizontal groups. Data on the horizontal group load dis- 
tributions during individual busy hours are not available for the test 
discussed in Section V and would be impractical to obtain on an on- 
going basis in any office. In the absence of data from which one could 
deduce directly the blocked-call peakedness, a treatment has been 
made using some results which R. I. Wilkinson obtained in the course 
of formulating his "equivalent random" method of characterizing 
overflow traffic. 11 

Wilkinson 11 assumes that traffic arrives and departs according to a 
birth-and-death process in which the arrival rate is increased when- 
ever the number of calls in the system exceeds a nominal number and is 
decreased whenever the number of calls in the system is less than this 
nominal number. The equilibrium distribution of the number of calls 
in the system resulting from these assumptions is shown to be negative 
binomial and, hence, this distribution is completely determined by its 
mean and variance (or by its mean and peakedness). Wilkinson then 
shows that negative binomial distributions rather closely approximate 
true overflow distributions under a number of different conditions. 

The effects of peakedness of the blocked-call stream on dial-tone 
delay were explored using the limiting analytical model of Section II 
by replacing the (truncated) Poisson distribution of blocked calls by 
a (truncated) negative binomial distribution of blocked calls. (Note 
that this does not require recomputing the conditional delay distribu- 
tions.) In this manner, dial-tone delay distributions were calculated 
and compared with the observed dial-tone-delay distributions dis- 
cussed in Section V, assuming peakedness values of 2 and 4. The re- 
sulting delay curves approximated the higher delays observed for 
given occupancy parameters. These calculations indicate that much 
of the observed variability in dial-tone delays under nearly identical 
load can be explained by peakedness of the blocked-call load. 

It is not difficult to modify the birth-and-death model discussed in 
Appendix B so as to accommodate peaked blocked-call input using 
Wilkinson's method of representing this input. 12 

Thus, when a negative binomial blocked-call distribution is used in 
the limiting model, the resulting distribution of the good-call dial- 
tone delay may be regarded as the limit of a sequence of good-call- 
delay distributions corresponding to models with blocked-call input 
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having the peaked form suggested by Wilkinson. 11 The limit is taken 
as the bad-call-arrival rates approach zero and the mean bad-call- 
waiting times approach infinity with their products approaching posi- 
tive constants. Incorporation of the negative binomial distribution in 
the limiting model is accomplished simply by replacing the terms of 
the truncated Poisson distribution in eq. (13) with the corresponding 
terms of the negative binomial distribution, 13 using any convenient 
truncation. 

V. COMPARISON OF CALCULATED AND OBSERVED DIAL-TONE 
DELAY DISTRIBUTIONS 

5.1 Summary 

Theoretical dial-tone-delay distributions, calculated using the ana- 
lytical model discussed in Section II, were compared with dial-tone- 
delay measurements made in a field test. The test was conducted in 
a No. 5 crossbar office with 60 llfs, 4 dtms, 68 dial-pulse (dp) ors, 
and 68 multifrequency (mf) ors. The data discussed in this section 
are from the time period February through April, 1974. 

The main conclusions of this study are that (i) the calculated delays 
generally fall in the midrange of applicable data, (ii) there is a large 
variability in observed dial-tone delays measured under nearly 
identical levels of dtm occupancy, second-failures-to-match, and 
percent of all-ORs-busy, and {Hi) the field data show a clear increase 
in the ratio of waste dtm usage to total dtm usage as frame load 
increases. The observed amount of increase in waste dtm usage agrees 
with theoretical predictions. 

In this section, the manner in which the dial-tone-delay measure- 
ments were taken is discussed and the method used to obtain the 
calculated delays is outlined. Next, some of the sources of variability 
in No. 5 crossbar dial-tone delays are identified and an explanation is 
given as to why a large variability in observed delays should be ex- 
pected. Two data plots are given indicating, respectively, the effects 
of frame load on waste dtm usage due to second-failures-to-match and 
the effects of frame load on incoming-first-failures-to-match. Finally, 
the conclusions of the study are discussed. 

5.2 Methods of obtaining the calculated and observed delay distributions 

Hourly dial-tone-delay measurements were made by placing ap- 
proximately 900 test calls per hour, using a standard 3-s dial-tone- 
delay testing machine which had been modified to record the pro- 
portion of test calls with delays exceeding X seconds, for X = 0.5, 
1.0, 1.5, 2.0, 2.5, and 3.0. 
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The observed actual dtm occupancy p b for a single hour is expressed 
in terms of the measured parameters by the formula : 

A = (Total dtm peg count) (0.015) + (Measured seconds of dtm usage) 
Pb (4) (3600) 

(16) 

The term in eq. (16) involving the dtm peg count is to adjust for the 
seconds of dtm usage which dtm usage measuring devices do not 
record. 

The (adjusted) dtm holding time for each hour was computed by 
A. R. Thorne from the dtm peg count, the all-ORs-busy peg count, 
the adjusted dtm usage [the latter of which comprises the numerator 
of eq. (16)], and an additional "light-traffic adjustment" (used when- 
ever the observed actual dtm occupancy is below 0.80). Most of the 
adjusted dtm holding times are between 0.28 s and 0.31 s. A dtm 
holding time of 0.30 s is assumed in the theoretical dial-tone-delay 
distributions discussed in this section. 

To compare predicted dial-tone delays with measured dial-tone 
delays, it was necessary to infer the observed increment in dtm occu- 
pancy due to dtm second-failures-to-match (dtm2ftm). This incre- 
ment, denoted by A, is taken to be 



. (Total dtm2ftm peg count) (HBC + 0.015) 

(4) (3600) ' (17) 



where 



HBC = dtm holding time during a second-failure-to- 
match (seconds). 

A value of 0.40 s is used for HBC, based on data obtained during an 
earlier dial-tone-delay field test. 6 The observed good-call dtm occupancy 
p u is defined to be 

Pu = h - A. (18) 

Dial-tone-delay distributions were calculated, using the analytical 
model, for a range of values of actual dtm occupancy p b and good-call 
dtm occupancy p u , where the parameters p b and p u are defined in 
Section II. The specific manner in which the distributions were calcu- 
lated is discussed below. The results of these calculations were tabu- 
lated into a set of dial-tone-delay distributions, indexed by the pairs 
(pb, p u )- To compare the calculated and observed dial-tone delays, 
measured dial-tone-delay distributions from data-collection hours with 
similar values of p b and p u were plotted on a graph along with one 
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theoretical dial-tone-delay distribution selected from the tabulation 
discussed above, such that pb tt p& and p u ~ Pu hold for the values 
/3 6 and p u corresponding to the observed delays shown on the graph. 

These graphs are shown in Figs. 3, 4, and 5. Table I lists data per- 
taining to each of the observed delay curves on each of the graphs. 
At the top of each graph are listed the actual dtm occupancy pb and 
the good-call dtm occupancy p u used for the theoretical dial-tone-delay 
curve (the solid line) on the graph. Also listed are the ranges of the 
Pb and p u values corresponding to the observed dial-tone-delay curves 
on the graph. The plotting symbols on the graphs indicate measured 
dial-tone delays. The dotted lines are smoothing curves fitted to the 
measured delays by the computer plotting routine used to draw the 
graphs. 

The predicted dial-tone-delay distributions were obtained in several 
steps. First, values of the good-call origination rate per llf Xi were 




2 3 4 

SECONDS OF DIAL-TONE DELAY 



Fig. 3 — Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.75 (data: 0.740 to 0.758). Good-call occupancy in calculations = 0.70 
(data: 0.684 to 0.697). 
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Fig. 4— Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.82 (data: 0.795 to 0.837). Good-call occupancy in calculations = 0.76 
(data: 0.744 to 0.782). 



computed so as to produce good-call dtm occupancies of p u = 0.60, 
0.66, 0.68, 0.70, 0.72, 0.74, 0.76, 0.78, and 0.80. The values of Xi were 
obtained numerically from the formula p u = \ t NT/2p given in Sec- 
tion III. Note that /;, the erlang C probability of all-ORs-busy, is a 
function of Xi, whereas .V and T are constants. For the office in which 
the test was conducted, X = ,30 and, as discussed above, T = 0.30. 
In computing p, an average or holding time of 10.25 s was assumed 
based on data from the test. In all cases, the calculated values of p 
were greater than 0.99, so the all-ORs-busy condition has a calculated 
probability of less than 0.01 under the conditions to which these 
distributions apply. (The observed fractions of all-ORs-busy were 
also below 0.01 during most of the hours of the test. Hence, ors do 
not appear to have caused much of the dial-tone delay observed in 
the test.) 

Next, the conditional delay distributions corresponding to these 
parameters were computed using the analytical model. Using eq. (12), 
values of x (total erlangs of blocked calls) corresponding to a range 
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0.006 



12 3 4 

SECONDS OF DIAL-TONE DELAY 



Fig. 5— Calculated and observed dial-tone delays. Actual occupancy in calcula- 
tions = 0.82 (data: 0.795 to 0.837). Good-call occupancy in calculations = 0.76 
(data: 0.744 to 0.782). (Note that a blocked-call peakedness of 4 is assumed in the 
calculations shown above. Figures 3 and 4 are based on an assumed blocked-call 
peakedness of 1.) 



of values of p b were computed for each fixed value of p u . For each such 
pair of p b and p u , the good-call dial-tone-delay distributions were then 
computed by (13) using the value of x, obtained as discussed above, 
and the conditional delay distributions corresponding to p u . The re- 
sult of these calculations is a set of good-call dial-tone-delay distribu- 
tions indexed by the pairs (p&, p u ). 

These good-call dial-tone-delay distributions include the indirect 
effect of bad calls in that they reflect the increased dtm congestion 
produced by the bad calls. As discussed earlier, the direct effect of bad 
calls is expected to be small and therefore may be accounted for in a 
somewhat approximate manner. To reduce the number of variables 
that need to be considered, a single blocking probability is used in 
lieu of integrating the blocking probability formula over an assumed 
distribution of expected horizontal group loads. For the purpose of 
computing an average blocking probability pb, the average waiting- 
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time-until-defection of a bad call (bcwt) is assumed to be 10 s.* A 
consequence of this assumption is that the total bad-call origination 
rate corresponding to x total erlangs of bad calls on two dtms is x/10 
bad calls per second. The total origination rate of good calls on all 
30 llfs is 30Xi. The average blocking probability pb corresponding to 
the occupancy pair (pj, p„) is then defined to be the ratio of the bad- 
call origination rate to the total origination rate of good and bad calls. 
Thus, 

PB = ,/(10)+30X, ' (19) 

where x and Xi are the values corresponding to the occupancy pair 

(Pb, Pu)- 

The conditional delay distribution of bad calls should be nearly 
exponential with mean delay ht/(10 — a), where ht is the average 
call holding time and a is the carried horizontal group load in erlangs. 1 
For each occupancy pair (p&, p u ), the theoretical dial-tone-delay dis- 
tribution for all calls is then given by 



P(D > t) = (1 - pb)P c (D > + pb jexp ' 



exp I - 



HT 



, (20) 



where P G (D > t) is the good-call dial-tone-delay distribution cor- 
responding to the occupancy pair (pb, p u ). 

In the theoretical dial-tone-delay curves shown in Figs. 3, 4, and 5, 
a = is used in eq. (20). The effect on P(D > 3 s) of setting a = 0, 
rather than using a more nearly correct value for each graph, is less 



'The value of bcwt = 10 s used in calculating pb is consistent with the mean- 
time-to-defection H b = 30 s used in the horizontal group blocking calculations in 
Section III. The difference between the numerical values arises because these two 
parameters are defined differently. For the purpose of calculating pb, it is assumed 
that a bad call arrives, remains waiting for dial tone an exponential length of time 
with mean bcwt, and then defects. In the horizontal group blocking model discussed 
in Section III, a call which finds all 10 line links busy may either defect or may 
eventually obtain an idle line link. The call contributes to the bad-call load during 
the time that it remains waiting for one of the line links to become available. Thus, 
the quantity in the horizontal group blocking model corresponding most nearly to 
the parameter bcwt is the mean time until a call that finds all 10 line links busy 
either obtains a line link or defects. Based on delay calculations for the horizontal 
group blocking model discussed in Section III, and assuming a mean-time-to-defec- 
tion of 30 s for calls that do not obtain a line link and a mean line-link holding time 
of 150 s, the mean time until a call that finds all 10 line links busy either obtains 
a line link or else defects is calculated to be between about 10 and 12 s for horizontal 
group carried loads in the range of interest. Thus, it is reasonable to take bcwt = 10 s. 

* This expression for the conditional delay distribution of bad calls should over- 
estimate their delays somewhat because the expression does not account for finite- 
source effects. Since the fraction of calls experiencing these delays (i.e., the fraction 
of bad calls) is typically less than 0.01, the effect on 3.0-s dial-tone-delay probabilities 
of neglecting finite-source effects is typically less than 5 X 10~ 4 . For this reason, 
these effects are neglected in eq. (20). 
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than 8 X 10~ 4 in all cases. This error is in the opposite direction from 
the error, of comparable magnitude, resulting from not including 
finite-source effects in the bad-call delay distribution. 

Since the analytical model is a discrete time model with a time step 
of T seconds, P(D > kT) = P[D ^ (k + l)r] for k = 1, 2, • • -. In 
Figs. 3, 4, and 5, T = 0.30. Because of the way that dial-tone-delay 
measurements are taken, a call which is recorded as having a delay 
greater than t seconds may actually receive dial tone within a few 
milliseconds after time t. Thus, in comparing the theoretical dial-tone 
delays with the observed dial-tone delays, the observed fraction of 
dial-tone delays greater than t seconds is taken as representing the 
observed fraction of dial-tone delays greater than or equal to t seconds. 
The theoretical delay curve plotted is the curve P(D ^ kT) for 
k — 1, 2, • • •, interpolated so as to produce a smooth curve. 

5.3 Sources of variation in observed No. 5 crossbar dial-tone delays 

Dial-tone delay in No. 5 crossbar is influenced by a number of 
factors capable of producing a large variation in delays measured in 
different hours within the same office under very similar conditions 
of dtm occupancy, percent all-ORs-busy, and second-failures-to-match. 
As discussed in Section IV, much of the variability in dial-tone delay 
measured under very similar load conditions can be explained by 
differences in peakedness of the blocked call stream. Whenever most 
of the blocked-call load comes from a small number of extremely over- 
loaded horizontal groups, the blocked-call stream should have a peaked- 
ness greater than one. When a large number of moderately overloaded 
groups contribute to the blocked-call load, the blocked-call stream 
should be approximately Poisson (peakedness equal to one). Thus, 
differences in the individual busy-hour-load balance would be expected 
to produce different amounts of blocked-call peakedness, which in 
turn can account for appreciable differences in dial-tone delay mea- 
sured under nearly identical average load conditions. For example, 
Fig. 5 shows dial-tone delays calculated assuming a blocked-call peaked- 
ness of 4 for the same conditions shown in Fig. 4, which is based on a 
blocked-call peakedness of 1.* The calculated delay curve in Fig. 5 fits 
the top two observed delay distributions rather closely. 

Some additional identifiable sources of dial-tone-delay variation 
under similar load conditions are: (i) within-hour trends in traffic, 
(ii) nonstandard (and possibly erratic) gating caused by improper 
functioning of the master traffic controller circuitry, (in) dtm prefer- 

* The method by which blocked-call peakedness is treated in the model is dis- 
cussed in Section IV. 
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ence for calls from a small subset of lines on each horizontal group, 
(iv) variation in dtm first-failures-to-match, and (v) competition be- 
tween dtms and completing markers for line-link connectors. 

The first of these sources should produce effects similar to those of 
blocked-call peakedness. The second source may cause nonuniform 
congestion. The third source is predicted to result in a slight outward 
shift in the delay curve. The fourth and fifth sources should be reflected 
in increased measured dtm holding time and in increased dtm occu- 
pancy. Although approximate allowances can be made for the average 
congestion increase produced by some of these phenomena, no quanti- 
tative estimate is available for their total contribution to hourly 
variation in dial-tone delay. 

In addition to identifiable sources of dial-tone-delay variation, simu- 
lation studies indicate that there can be an appreciable residual varia- 
tion in simulated hourly No. 5 crossbar dial-tone delays obtained in 
different runs with identical input parameters (and, hence, with identi- 
cal expected load conditions).* 

Figure 6 shows four dial-tone-delay distributions obtained using 
the gating simulation model. In this model, the blocked-call stream 
is Poisson. These distributions were produced by simulating four 
individual hours, using identical input parameters. The delay distribu- 
tions shown are for the calls that did not encounter horizontal group 
blocking and are based on the total number of such calls processed 
during the hour. The set of four 3-s dial-tone delays has a coefficient 
of variation of 0.28 and a mean of 0.069. (The coefficient of variation 
is the ratio of the standard deviation to the mean.) Plotted on the 
graph along with the delay curves are error bars indicating the 2-sigma 
limits of 0.034 and 0.107 associated with the above mean and coeffi- 
cient of variation. 

Actual dial-tone-delay measurements are based on test calls. During 
a given busy hour in a typical No. 5 crossbar office, approximately 
900 test calls are made on a fixed set of 60 (out of 600) horizontal 
groups. The use of test calls introduces sampling error, which is not 
represented in the distributions shown in Fig. 6 and which would 
have an associated coefficient of variation of about 0.12 for the parame- 
ters applicable to Fig. 6. 

5.4 Data on frame load effects 

Figure 7 is a data plot of line-link frame load versus waste dtm usage 
due to second-failures-to-match based on data from the test. The 
quantity "dtm usage fraction due to 2ftms" shown on the ordinate is 



* This result is one of the main conclusions of an earlier No. 5 crossbar dial-tone- 
delay simulation study conducted by S. Halfin. 3 
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Fig. 6 — Simulated No. 5 crossbar dial-tone delays. Distributions are based on simu- 
lation of four individual hours with identical inputs. 



defined by 
dtm usage fraction due to 2ftms = Total dtm usage due to 2ftms 



Total dtm usage 



A 

pb 



Figure 8 is a data plot of line-link frame load versus incoming-first- 
failure-to-match (iffm) based on data from the test. In each figure, 
each of the data points represents data from 1 h. The average actual 
dtm occupancy (averaged over all the data-collection hours with 
measured line-link frame loads of 1100 ccs/llf or more) is 0.54; the 
dtm occupancy range is from less than 0.40 to 0.84. 

Existing theory indicates that, for a given office configuration (in- 
cluding a given junctor pattern), iffm is directly and primarily de- 
pendent on frame load. This conclusion is borne out by Fig. 8, which 
demonstrates a well-defined trend (with a moderate amount of data 
scatter) of increasing iffm with increasing frame load. 

The extent to which frame load affects dtm usage and dial-tone 
delay in any given hour depends, in an indirect way, on several differ- 
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Fig. 7 — Effect of frame load on increased dtm usage due to 2ftMb. The assumed 
dtm holding time during a 2ftm = 0.40 s. 



ent variables, including dtm occupancy and the distribution of carried 
load among horizontal groups. Because of this dependence, the nature 
of which has been discussed more fully in previous sections, the data 
scatter in the plot of waste dtm usage fraction versus frame load (Fig. 
7) is much larger than in the plot of iffm versus frame load (Fig. 8). 

The curve labeled "theoretical" in Fig. 7 was calculated using the 
limiting analytical model [eq. (12)] and the horizontal group blocking 
model discussed in Section III. In these calculations, the blocked-call 
load was assumed to be Poisson (peakedness equal to 1). 

To obtain the theoretical curve, it was first necessary to determine 
what calling rates should be assumed in the calculations. The calling 
rates were inferred from the data upon which Fig. 7 is based by first 
using linear regression (least squares) to express the observed good-call 
dtm occupancy, p u , as an empirical function of the observed frame load, 
ccs/llf. For each value of the frame load, the calling rate was taken 
to be the particular calling rate corresponding to the least squares 
value of p u , assuming an average dtm holding time of 0.30 s and a 
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Fig. 8 — Effect of frame load on incoming first-failure-to-match. 

zero probability of all-ORs-busy. The justification for the assumptions 
regarding the dtm holding time and the probability of all-ORs-busy 
is given earlier in this section. 

Next, the horizontal group blocking model discussed in Section III 
was used to compute the expected blocked-call offered load correspond- 
ing to each value of frame load and (empirically associated) calling 
rate. The office distribution of carried horizontal group load was 
assumed to be normal with group-to-group coefficient of variation 
inferred from the office horizontal group load distribution using the 
method discussed earlier. The calculated values of offered blocked- 
call load were then used in eq. (12) to compute the theoretically pre- 
dicted fraction of waste dtm use, (p b - p u )/p b , corresponding to these 
frame loads and calling rates. 

Since both the theoretical and the observed values of waste dtm 
use depend not only on frame load but also on the calling rate, neither 
the data plotted in Fig. 7 nor the theoretical curve shown on the figure 
should be regarded as being applicable to calling rates or load balance 
conditions other than those upon which this figure is based. 
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The purpose of Fig. 7 is to illustrate — for the conditions of calling 
rate, frame load, and load balance represented in this study — an 
empirical relationship between frame load and the ratio of waste dtm 
usage to total dtm usage with increasing frame load and to show that 
this empirical relationship is consistent with predictions of the theo- 
retical models. 

5.5 Conclusions 

The conclusions which the author has drawn from this comparison 
of theoretical and observed dial-tone-delay distributions are as follows : 

(i) Where two or more observed delay distributions appear on the 
same graph, the theoretical delay distributions usually fall 
approximately midway between the maximum and minimum 
observed delay distributions. This indicates that, in the 
(actual) dtm occupancy range spanned by these data (dtm 
occupancies up to about 0.84), the analytical model shows 
good agreement with the data. 
(it) Much of the large variability in observed 3-s dial-tone delays 
measured under nearly equal conditions of dtm occupancy, 
second-failures-to-match, and percent all-ORs-busy can be ex- 
plained by assuming the blocked-call stream to have different 
peakedness values (ranging from 1 to about 4) in different 
hours. High blocked-call peakedness, illustrated by Fig. 5 in 
which the peakedness is taken to be four, would be expected 
whenever most of the blocked-call load comes from a small 
number of extremely overloaded horizontal groups. Low 
blocked-call peakedness, illustrated by Fig. 4 in which the 
peakedness is taken to be 1, would be expected whenever a 
large number of moderately overloaded horizontal groups con- 
tribute more or less equally to the blocked-call load. 
(Hi) Simulation results indicate that there should be a large residual 
variability in 3-s dial-tone delays measured in different 
(simulated) hours under identical expected load conditions. 
This variability, illustrated in Fig. 6, is in addition to the 
variability due to blocked-call peakedness discussed above. 
A third nonnegligible source of variability in observed dial- 
tone delays is the sampling variability associated with the use 
of test calls to measure delays. 
(iv) For the conditions of calling rate, frame load, and load balance 
represented in this study, the observed increase in the ratio of 
waste dtm usage to total dtm usage as frame load increases is 
consistent with theoretical predictions. This is illustrated in 
Fig. 7 and discussed in Section 5.4. 
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APPENDIX A 

Computation of the Two-Marker Stationary Distribution for k^1 

The purpose of this appendix is to prove that the limit r defined 
by (10) exists and that the quantities Vj = <xj + rfa, satisfy (21), 
where a ; and 0,- are defined by (9). 

Vj ^ 0, all j ^ 

v = 1; 

Vj = E vf^i, j), all j fc 

i=0 

£«/<•■ (21) 

i=o 

We first show that there exists exactly one positive number r such 
that Vj = a, + r@j satisfies (21). Then we show that the limit in (10) 
exists and is equal to this number. 

When \ < 2p, we know that the stationary distribution P k exists 
and P k (j) > for each j. Let irj = P k (j)/P k (0). Then for r = P k (l)/ 
P k (0), we have 

no = 1, 
in = r, 

*i = E nPk(i, j), all j ^ 0. 

i=0 

Now let Vj = ctj + r/3y, with r as defined above. Then, it follows from 
(9) that 

v = 1, 

vi = r, 

^■=E^,i), allJ^O. (22) 

i=0 

Since v — tr and v x = n and since the vectors v and n both satisfy 
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the same recurrence relation, which has the property that the jth. 
term for each j ' ^ 2 is uniquely determined by the Oth and 1st terms, 
we see that v, = t } - > for all j ^ 0. Since the Markov chain is posi- 
tive recurrent for X/2p < 1, it follows from (22) and Karlin 14 that 
L"-o«i <°°- Hence, v satisfies (21). Thus, there exists at least one 
value of r for which v 3 - = a, + rfr satisfies (21). If r' is any number 
such that v'j = ay + r% satisfies (21), then by uniqueness of the 
stationary distribution, r' = vM = P*(1)/P*(0), so r' = r. Thus, 
exactly one such r exists. 

Since a,- + rjSj > for each /, it follows that (10a) holds. Hence, 
the increasing sequence ra n is bounded above by r and the decreasing 
sequence M„ is bounded below by r. So both m = lim„ ra„ and 
M = lim„ M„ exist and satisfy 

m << r g M, (23) 

^n ^ m, 4) 

M £M n . 

If either of the inequalities in (23) were strict, then m < M and, in 
view of (24), each x in the interval m < x < M would satisfy 

for each n, from which it follows that Vj = a } - + x0 3 - satisfies (21). 
This contradicts the uniqueness of r and proves that m = M = r. 

It is easy to extend the above result so as to give a method for 
computing the stationary distribution of any irreducible positive re- 
current Markov chain on the nonnegative integers such that a fixed 
positive integer n exists for which P(j + n , j) > and P(j + n, j) 
= for all j and all n > n . 

APPENDIX B 

Convergence in Distribution to the Limiting Model 

In the limiting model, the number of bad (i.e., blocked) calls in the 
system is a time-independent, truncated-Poisson-distributed random 
variable. We will show that the steady-state distributions of good- and 
bad-call queue lengths and the steady-state good-call dial-tone-delay 
distribution in the limiting model are the limits of the corresponding 
distributions for a sequence of models in which the bad call queue 
lengths form birth-and-death processes. 

In the nth model, the bad-call arrival rate, denoted by X (n) , ap- 
proaches zero and the bad-call mean waiting time until defection, de- 
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noted by H< n \ approaches infinity in such a way that their product 
(total expected erlangs of bad calls) approaches a finite positive con- 
stant x. 

The physical motivation for considering the limiting model is that 
the actual situation is one in which blocked calls appear infrequently 
(relative to total call arrival rates) and tend to remain in the system 
for a long time (relative to dtm holding times). The mathematical 
motivation is that the limiting model is much easier to analyze than 
a model in which the bad-call arrival and departure processes are 
represented explicitly. 

Preliminary computations using a single marker version of the 
analytical model, which represented the bad-call arrival and departure 
processes explicitly, showed that different bad-call arrival rates and 
waiting times had little effect on dial-tone-delay distributions as long 
as the product of the bad-call arrival rate and bad-call waiting time 
remained constant. In addition, these distributions were all quite 
close to those obtained from the limiting form of the model. These 
results make sense intuitively because bad calls simply cycle through 
the system, absorbing some dtm uses while present, and eventually 
defect; thus, what should matter is mainly the distribution of the 
number of bad calls in the system at any time. In the version of the 
analytical model discussed in this appendix, the blocked call queue 
length distribution is shown to be a truncated Poisson with mean 
equal to the total erlangs of bad calls. 

Let K be the maximum possible number of bad calls that can be 
in the system at any time. In all computations using realistic No. 5 
crossbar busy-hour input parameters, the expected erlangs of bad 
calls have been low enough that the Poisson probability of more than 
five blocked calls being present in a two-marker system has been less 
than 10~ 3 . Hence, K may be taken to be 5. (Note that in an actual 
No. 5 crossbar office, the maximum number of bad calls that can be 
in the system at any time is trivially bounded above by the number of 
subscriber line terminations in the office.) 

Let S„(p) denote the saturation load of the nth model. Since no 
more than K blocked calls can be in the queue at any time, in the 
nth model, we have 2p = Sg*(p) ^ S n (p) ^ Sg>(p) = 2p for all n; 
hence, S„(p) = 2p. Thus, the saturation load of the nth model is the 
same as the saturation load of the limiting model. For the nth model, 
let 

X™ = number of good calls in the queue at time mT, 
Y (n) (mT) = number of bad calls in the queue at time mT, 
p n {k,l) = P{Y^l(m + l)r] = l\Y™(mT) = k). 
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In the nth model, Y M (t), the bad-call population size at time t, is a 
finite-state space birth-and-death process with birth-and-death rates 

X™ = X (n) ^ m ^ K - 1 
= otherwise 

u (n) _ m < m < K 

= otherwise. 

Since the queuing model is in discrete time, only the values of F (n) (i) 
at i = mT are of interest. The transition matrix p n (k, I) is given by 

p _ e TA n ^ 

where A„ is the infinitesimal generator matrix of the birth-and-death 
process F (n) . (SeeRef. 15.) Since lim,-* X< n > = Oandlim„H.« # (n) = «, 
it follows that lim,,-* A„ = 0. Hence, 

lim p»(fc, I) = 8 kt (25) 

n-»°o 

for each k, I = 0, 1, • • •, K, where 5*1 = for k ^ I and 8 kk = 1. 
We also have 

where Pjt(i, i) is given by (1) and (2). Since the nth model has satura- 
tion load 2p, we know that a stationary queue length distribution 
tnij, I) exists for X/2p < 1 and satisfies the equation 

ir»(i, I) - £ E *n(«, k)Vn(k, l)P k (i, j). (26) 

fc=0 i=0 

Let 7T„( • , I) be the marginal equilibrium distribution of the number of 
bad calls in the system. Then 



T»(-,0 = £Tn(y,Z) (27) 

j-0 



and by (26) 



».(■, = EE ».(*, fc)Pn(fc, f P*(i, j) 

A;=0 t=0 /— 

= !>„(-, fc)p»(M)- (28) 

fc=0 

Hence, ir„( • , /c) is the stationary distribution of a Markov chain whose 
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transition matrix is p n (ft, I). Using standard theorems on birth-and- 
death processes (Ref . 10) it follows that, for all values of I = 0, 1, • • • , K, 

Tn(-,k) =limP[F<«>(0 = ft|7C">(0) = Z] 



«-»« 









A;! ~~ 








= otherwise, 


where 






X n = X (n) ff< B > 


and 






K ~t 

C*nK = 2- Ti' 

k=0Kl 


We 


will 


now 


show that 

lim irJi. K) - CxKxk F 



(29) 



(30) 

for all ft - 0, 1, • • •, K and all j = 0, 1, • ■ ■ , where P k is the stationary 
distribution of good-call queue length in the model with ft bad calls 
permanently present in the system. The right-hand side of (30) is the 
queue-length distribution for good and bad calls in the limiting model. 
Let X (n) denote the number of good calls in the queue for the nth 
model in equilibrium. Since the number of bad calls present in the 
/tth model is always less than or equal to K for all n, it is clear on 
intuitive grounds and can be proved rigorously using stochastic order- 
ing that 

P[X™ ^ jl*t Pk® (31) 

for each j. 

To prove (30), it suffices to show that if t„> is any subsequence of 
7r„ for which lim„- w»*(j t ft) ■ ir(j, ft) exists for each j and ft, then 



and 



(*) T(j,k) = Zir(i,k)P k (i,j) 



-Art 



CxkX 



(**) £7T(j,ft) = ^- 

j=o ft! 



To see that (i) and (ii) are sufficient, note that, by uniqueness of 
the stationary distribution, P k , (i) and (ii) imply 

*U,® = Cj ^P k (J)- (32) 
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Since ^ ir n ( j, k) ^ 1 for all j and k, every subsequence t„> contains 
a further subsequence ir n " for which the limit lim„" ic n "(j, k) exists 
for all j and k. In view of the above discussion, all of these subse- 
quences have the same limit, namely the right-hand side of (32). 
Equation (30) follows. 

We now prove (i) and (ii). Condition (i) follows immediately from 
(25) and (26). To see that (ii) holds, we proceed in two steps. First, 
note that by (27), (29), and Fatou's lemma, 

for each A;. If any of the above inequalities were strict, then 

e i *a fc) < i. 

*=o y=o 
Hence, to prove (ii) it suffices to show 

t t *(h *) ^ 1, 

t— i=0 

i.e., we must show that, in the limit, no probability mass escapes to 
infinity. Let e > and choose jo such that 

E PkU) < 6. 
3=30 

Then, using eq. (31), 

f E *(*, k)^T. Y, »(f, *) 

i=0 jfe=o i-0 A=0 

io JC 

= lim E 2Z T»"(»i fc) 

n" t=0 k=0 

= lim P[X<»"> g >] 

n" 

^ E P*(0 

^ 1 - e. 

Thus, (ii) holds and eq. (30) follows. 

It can be readily shown, using the above results together with 
standard theorems, that the delay probabilities in the nth model 
converge to those in the limiting model. (See Refs. 16 and 17.) 

1792 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1975 



REFERENCES 

1. S. Halfin, unpublished work. 

2. F. R. Wallace, personal communication. 

3. S. Halfin, "An Approximate Method for Calculating Delays for a Family of 

Cyclic Type Queues," B.S.T.J., this issue, pp. 1733-1754. 

4. H. Kushner, Introduction to Stochastic Control, New York: Holt, Rinehart, and 

Winston, 1971, p. 211. 

5. W. S. Hayward, Jr., "Traffic Engineering and Administration of Line Concen- 

trators," Congress Record, Paper No. 23, Second International Teletraffic 
Congress, The Hague, 1958. 

6. J. G. Kappel, personal communication. 

7. C. Clos and R. I. Wilkinson, "Dialing Habits of Telephone Customers," B.S.T.J., 

31, No. 1 (January 1952), p. 46, last paragraph. 

8. J. Riordan, Stochastic Service Systems, New York: John Wiley, 1962, pp. 109-112. 

9. Ref. 7, pp. 32-67. 

10. S. Karhn, A First Course in Stochastic Processes, New York: Academic Press, 

1968, p. 194. 

11. R. I. Wilkinson, "Theories for Toll Traffic Engineering in the USA," B.S.T.J., 

35, No. 2 (March 1956), pp. 421-514. 

12. Ref. 11, p. 453. 

13. Ref. 11, p. 454, eq. (19). 

14. Ref. 10, theorem 3.2. 

15. Ref. 10, pp. 206-208. 

16. P. Billingslej', Convergence of Probability Measures, New York: John Wiley, 

1968, p. 224. 

17. H. Royden, Real Analysis, 2nd Edition, New York: MacMillan, 1968, p. 232. 



DIAL-TONE DELAY 1793 



